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随 着 表面 科学 的 发 展 , 人 们 需要 更 深入 了 解 材 五 
料 表 面 的 物性 . 材料 表面 的 表面 能 和 功 函数 @ 是 。 帮助 认识 材料 表面 的 稳定 性 、 表 面 化 学 反应 发 生 的 
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摘 要 对 Be, Mg, Ti 和 Zr4 种 六 方 结构 的 纯 金属 进行 不 同 表面 的 功 函 数 和 表面 能 的 计算 . 结果 表明 , 对 于 这 4 种 金属 而 言 ， 
不 同 表面 取向 之 间 的 功 函数 和 表面 能 差别 很 大 , 其 中 (0001) 表 面 的 功 函数 最 高 , 表面 能 最 低 , (0111) 表 面 的 功 函 数 最 低 , 表面 
能 最 高 ; 同族 的 2 种 金属 , 其 不 同 表面 的 功 函 数 @ 变 化 规律 相似 , 表现 为 : Q(0001D)>G(0221)>G(1010)>G@(1231)>G(1230)> 
@(0111), 而 表面 能 y 的 变化 规律 完全 相反 , 表现 为 : KX0001)<X0221)<X1010)<X1231)<X1230)<X0111); 同族 的 2 种 金属 , 原 
子 序数 大 的 金属 , 其 同一 表面 取向 的 功 函数 和 表面 能 都 要 相对 另 一 种 金属 的 小 . 
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ABSTRACT Surface energy and work function are important for understanding metal surface performance. 
There have been many reports concerning these problems for bcc and fcc metals. However, for hex metals, there 
are few reports, especially for some surfaces with higher indices. This work focuses on the calculation of surface 
energies and work functions for Be, Mg, Ti and Zr. Many different surfaces with higher indices are calculated, in- 
cluding (0001), (0221), (1010), (1231), (1230) and (0111). The results show that the surface energies and work 
functions have a strong relationship with each other. The work functions of these surfaces follow such sequence: 
D0001)> T0221)> B1010)> 1231)> (1230)> B0111), corresponding their surface energies: A0001)<X0221)< 
M1010)<X1231)<X1230)<X0111). What's more, the results also show that, for these metals in the same group, the 


small ones in atomic number such as Be and Ti, compared with Mg and Zr, have larger values of work functions 
and surface energies. In addition, the (0001) surface has the largest value of work function but smallest value of sur- 
face energy, which is in good agreement with experiment outcomes. 


KEY WORDS work function, surface energy, first-principle 


究 其 表面 性 质 十 分 基础 且 重要 的 数据 . 它们 可 以 
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点 基础 研究 发 展 计划 资助 项 目 2011CB606403 等 众多 理化 特征 . 因此 , 对 于 材料 表面 能 和 功 函数 
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的 研究 十 分 重要 . 
Skriver 等 中 和 Andersen 等 2 利用 TB-LMTO 
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(tight-binding linear muffintin orbital) 方 法 , 系统 研究 1 计算 方法 
了 材料 的 表面 能 和 功 函 数 , 但 他 们 考虑 的 多 为 材料 1.1 表面 能 


的 常规 表面 , 对 于 同一 金属 高 指数 的 表面 功 函数 和 
表面 能 未 做 过 系统 研究 . Rohwerder 和 Turcu5 通 过 
扫描 kelvin 探 针 力 显微镜 (SKPFM) 直 接 测 量 材料 表 
面 的 电势 数据 , 但 由 于 方法 本 身 的 限制 , 待 测量 节 
表面 很 难 达 到 理想 的 状态 , 最 终 得 到 的 只 是 一 个 综 
合 表面 各 因素 的 结果 , 极 难得 到 纯净 表面 的 准确 的 
相关 数值 , 据 此 难以 系统 研究 金属 不 同 表 面 的 表 再 
能 和 功 函 数 . 
基于 密度 泛 函 理论 研究 材料 表面 功 函 数 和 表 
面 能 的 工作 有 很 多 . 利用 密度 泛 函 理论 计算 方法 ， 
能 够 在 合理 的 简化 表面 模型 后 , 在 短 时 间 内 进行 大 
量 的 不 同 表 面 类 型 的 计算 , 甚至 也 可 以 研究 一 些 实 
际 中 很 难 被 观察 到 但 理论 上 可 以 存在 的 表面 . 
Kokko 等 9 系统 研究 了 Li 的 (100), (110) 和 (111) 3 个 
表面 的 功 函数 和 表面 能 , Singh-Miller 和 Marzari"" 系 
统 研究 了 Al, Pd, Pt, Ti 等 金属 的 几 个 常见 低 指 数 表 
面 (001), (110), (L1D 和 (0001) 的 功 函 数 和 表面 能 ; 
Skriver 和 Rosengaardm 研 究 了 一 系列 纯 金 属 的 密 排 
面 的 表面 能 和 功 函 数 , 认为 金属 密 排 面 的 的 功 函 
数 和 表面 能 与 其 原子 序数 有 关 , 而 且 对 于 fcc 金 
属 存在 如 下 关系 : @(111)>@D(100)>@D(110), X110)> 
MA100)>X111). 

上 述 研究 都 是 基于 金属 的 常见 表面 , 如 立方 晶 
系 的 (001), (011), (111) 表 面 以 及 六 方 品系 的 (0001) 
表面 . 由 于 六 方 品系 晶 胞 的 ww2 轴 唱 格 参数 与 c 轴 唱 
格 参数 差别 很 大 , 导致 不 同 晶 面 指数 的 表面 理想 结 
构 较 为 复杂 , 极 大 增加 了 其 研究 难度 . 目前 还 很 少 
见 到 对 于 六 方 晶 系 不 同 晶 面 指数 表面 的 研究 报道 ， 
大 多 数 是 关于 其 (0001) 晶 面 的 相关 研究 工作 p 

近期 Wang 等 对 于 bcc 和 fee 体系 若 干 纯 金属 
的 一 系列 不 同 表 面 的 表面 能 和 功 函 数 进 行 了 研究 ， 
发 现 对 于 bcc 金 属 (Nb, Mo, Ta 和 W), 功 函 数 遵 从 如 
下 规律 : @(110)>@(133)>@(120)>@(111)>@DG311)> 
(100); 对 于 fcc 类 型 的 金属 , 功 函 数 遵从 如 下 规 
律 : Q(11D)>G(100>G(C11D)>G123)>G(310)>G(C110). 
Tang 等 I 利用 第 一 原理 方法 对 六 方 晶 系 的 纯 Mg 若 
干 表面 的 表面 能 进行 了 研究 . 但 到 目前 为 止 , 关于 
其 它 一 些 六 方 晶 系 金 属 以 及 六 方 晶 系 金属 的 一 些 
其 它 晶 面 指数 的 相关 计算 研究 尚未 见报 道 . 


表面 能 可 以 很 好 地 反映 材料 表面 的 稳定 性 等 
言 息 , 目前 , 利用 计算 模拟 来 计算 表面 能 的 方法 有 
很 多 , 其 本 质 是 计算 构建 的 表面 模型 与 其 块 体 模型 
之 间 的 能 量 差 . 具体 计算 时 , 当 模 型 的 上 下 2 个 表面 
为 同一 种 表面, 表面 能 的 计算 如 下 式 所 示 m 
y= 2 (Els ~ NEvar) (1) 
式 中 , EY, 是 表面 模型 的 总 能 ，N 是 表面 模型 中 的 
原子 个 数 ，Eii 是 材料 块 体 中 单个 原子 的 能 量 , 8 是 
模型 表面 面积 . 在 建立 模型 时 , 要 充分 考虑 了 上 下 2 
个 表面 的 具体 情况 , 保证 计算 时 2 个 表面 状态 相同 ， 
为 同一 种 表面 . 
1.2 功 函 数 
功 函数 定义 为 从 材料 内 部 移出 一 个 电子 到 材 
料 表 面 外 一 点 所 需要 的 最 小 能 量 , 即 电子 摆脱 金属 
束缚 必须 做 的 功 . 可 用 下 式 表示 : 
D= Em rom (2) 


ermi 


材料 的 Fermi 能 级 . 

实际 上 , 利用 第 一 原理 计算 时 , 很 难得 到 
和 ,的 绝对 数值 , 不 过 能 够 同时 得 到 两 者 的 相对 
值 , 直接 计算 其 差 值 , 得 到 的 功 函 数 信息 与 真实 值 
是 相同 的 , 可 用 下 式 表述 : 


计算 时 , 天，,; 是 相对 变化 的 , 但 BF',, 改变 ,EE',， 
也 会 发 生 相 应 的 变化 . 此 时 式 (3) 得 到 的 功 函 数 与 式 
(2) 得 到 的 结果 在 数值 上 是 一 致 的 . 图 1a 给 出 了 计算 
功 函 数 时 的 模型 简 图 , 图 lb 给 出 了 对 应 于 此 模型 的 
势能 分 布 曲线 . 通过 计算 可 以 得 到 忆 ，; 和 已 ，，， 
计算 两 者 的 差 值 , 即 可 得 到 该 表面 模型 的 功 函数 . 
2 计算 模型 与 细节 
2.1 计算 模型 

本 工作 研究 的 对 象 为 Be, Mg,Ti 和 Zr 4 种 六 方 
晶 系 的 纯 金 属 , 分 别 计算 了 其 (0001), (0111), (0221)， 
(1010), (1230) 和 (1231) 6 种 不 同 表面 构 型 的 表面 能 
和 功 函 数 . 计算 采用 周期 性 边界 条 件 , 分 别 建立 了 
每 种 金属 的 6 种 不 同 表 面 的 超 胞 模型 , 然后 进行 相 


本 工作 基于 对 Be, Mg, Ti 和 Zr 4 种 六 角 唱 系 金 


Kr 


属 的 6 种 表面 进行 计算 研究 , 来 探讨 其 表面 的 功 函 
数 与 表面 能 的 相关 规律 . 


关 优化 和 计算 . 超 胞 模型 采用 的 是 横向 为 1x1 构 型 ， 
纵向 为 一 定金 属 层 数目 的 表面 模型 . 模型 中 金属 原 
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工 功 函数 计算 简 图 


Fig.1 Schematic diagram of the calculation of work func- 


tion 
(a) model picture 
(b) electrostatic potential curve of the model (Enmr 一 
calculated Fermi energy, 已 (sun 一 calculated vac- 
uum energy，c direction is vertical to the model 


Surface) 


(a) 


图 2 计算 时 Mg 的 不 同 表面 的 模型 示意 图 
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子 层 的 厚度 选取 根据 表面 构 型 的 复杂 程度 而 有 所 
不 同 , 但 都 优先 保证 了 模型 上 下 2 个 表面 的 对 称 性 
相同 , 以 确保 上 下 2 个 表面 反映 的 为 同一 种 表面 . 图 
2 给 出 了 计算 时 采用 的 Mg 的 6 种 不 同 表 面 的 构 型 
示意 图 . 其 它 3 种 金属 的 表面 模型 与 其 很 类 似 . 
2.2 计算 细节 

计算 采用 VASP (viennaab initio simulation pack- 
age) 软 件 包 0, PAW (projector-augmented wave) 方 法 ， 
PBE (Perdew-Burke-Ermmzerhof) 广 义 梯度 近似 "*, 平 
面 波 截断 能 选取 为 400 eV, 真空 层 厚度 统一 取 为 
1.6nm. 开 点 网 格 根据 表面 模型 厚度 选取 , 真空 层 方 
向 由 于 较 厚 ,区 点 基数 统一 取 为 1, 但 模型 K 点 总 数 
不 低 于 40 个 . 计算 时 , 相关 收敛 标准 为 系统 总 能 收 
敛 精度 小 于 1x10" eV, 作用 在 每 个 原子 上 的 力 小 于 
0.002 eV/ (nm.atomy). 
有 具体 计算 前 对 4 种 金属 的 品格 常数 进行 了 优化 
计算 , 并 与 实验 数值 "进行 了 比较 , 结果 见 表 1. 可 以 
看 出 , 计算 得 到 的 4 种 金属 的 a 轴 与 < 轴 的 蝇 格 常数 以 
及 晶体 的 轴 比 与 实验 值 之 间 的 差别 很 小 , 计算 值 与 实 


Fig.2 Models used in this calculation for (0001) (a), (0111) (b), (0221) (c), (1010) (d), (1230) (e) and (1231) (f) of Mg 


表 1 晶 格 常数 的 验证 


Table 1 Comparisons of the calculated crystal lattices with experimental crystal lattice constants 


Metal Calculation result Experimental value'" 
41nm c/nm c/a a/nm c/nm c/a 
Be 0.227 0.356 ls57 0.229 0.358 由 3 
Mg 0.319 0.519 1.63 0.321 0.521 1.62 
Ti 0.293 0.465 1.39 0.295 0.468 139 
Zt 0.324 0.518 1.60 0.323 0.515 1;39 
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验 值 之 间 的 误差 不 超过 1%, 在 可 接受 的 范围 内 . 
3 Be, Mg, Ti 和 Zr 不 同 表面 的 计算 结果 与 讨论 
对 Be, Mg, Ti 和 Zr 4 种 金属 的 (0001) 表 面 的 计 
算 结 果 与 文献 报道 的 计算 结果 以 及 实验 数据 "3 进 
行 了 比较 , 具体 见 表 2. 关于 其 它 表 面 的 计算 结果 ， 
由 于 数据 报道 很 少 , 在 此 未 作 比 较 . 从 表 2 可 以 看 
出 , 本 工作 的 计算 结果 与 文献 报道 的 计算 结果 很 接 
近 , 且 计 算 结果 与 实验 数值 差别 不 大 , 可 以 认为 本 
工作 中 计算 得 到 的 表面 能 和 功 函 数 的 数据 是 可 靠 
的 , 可 以 用 来 进行 相关 数据 分 析 . 
3.1 Be 和 Mg 的 数据 结果 与 讨论 
对 Be 和 Mg 相关 表面 的 计算 进行 了 计算 验证 ， 
构建 模型 的 金属 原子 层 厚 度 由 12~36 层 不 等 , 计算 时 
根据 不 同 表面 的 复杂 情况 而 调整 . 这 个 厚度 达到 其 
至 超过 了 目前 已 查阅 的 文献 报道 "中 所 采用 的 原子 
层 厚度 , 他 们 大 都 采用 13 层 左右 的 金属 原子 层 厚度 
进行 建 模 . 本 工作 之 所 以 采用 比较 大 的 原子 层 数 建 
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达到 收敛 的 标准 . 不 同 厚 度 下 表面 能 收敛 得 很 好 ， 
eh 0.03 Ji 计算 得 到 的 功 函 数 收敛 也 很 好 , 
波动 小 于 0.05 eV. 实际 计算 时 , 由 于 选取 的 原子 层 
数 较 多 , 每 个 层 数 模 型 下 得 到 的 表面 能 和 功 函 数 都 
已 稳定 , 每 一 层 数 下 的 计算 结果 都 可 视 为 收敛 , 因 
此 功 函 数 和 表面 能 的 具体 取 值 是 计算 收敛 后 的 几 
个 不 同 厚度 的 结果 的 平均 值 . 
图 4 给 出 了 Be 和 Mg 不 同 表面 的 功 函 数 和 表面 
能 的 计算 结果 . 可 以 看 出 , Mg 和 Be 不 同 表面 之 间 的 
功 函 数 和 表面 能 差异 很 大 . 对 于 Mg 而 言 , (0001) 表 面 
的 功 函 数 最 高 , 为 3.71 eV, (0111) 表 面 的 功 函数 最 低 ， 
为 3.48 eV, 不 同 表面 功 函 数 差异 达 0.23 eV; (0111) 表 
面 的 表面 能 最 高 , 为 0.84 Jm2, ( 0001) 表 面 的 表面 能 
最 低 , 为 0.56 J/m’, 表面 能 差异 达 0.28 J/m”. 对 于 Be 
而 言 , 规律 很 相似 , ( 0001) 表 面 的 功 函 数 最 高 , 为 
5.23 eV, ( 0111) 表 面 的 功 函 数 最 低 , 为 4.22 eV, 不 同 
表面 功 函 数 差异 可 达 1.01 eV; (0111) 表 面 的 表面 能 


模 , 主要 是 考虑 到 某 些 表面 构 型 偏 复杂 , 层 间 距 偏 
小 , 为 了 保证 模型 中 金属 部 分 的 厚度 达到 1.0 nm 磊 
, 因此 增加 了 原子 层 厚度 . 计算 过 程 中 所 有 结果 
和 数据 分 析 都 是 在 计算 数据 已 达到 设 定 的 收敛 标 
准 后 而 进行 的 ， pa 可 以 看 出 , 对 于 不 
同 的 表面 取向 , 选取 的 金属 层 厚度 不 同 , 但 都 已 经 


与 


最 高 , 为 2.83 J/m”, ( 0001) 表 面 的 表面 能 最 低 , 为 
1.76 J/m’”, 差异 达 1.07 J/m*. 2 种 金属 不 同 表面 之 间 
的 功 函 数 排列 规律 近乎 相同 , 且 功 函数 由 高 到 
低 的 顺序 为 : QC0001D)>GXK0221)>G(1010)>G(1231)> 
GO(1230)>G0111D); 2 种 金属 不 同 表 面 的 表面 能 的 
高 低 顺 序 也 完全 相同 , 但 表面 能 的 高 低 顺 序 刚好 


表 2 (0001) 表 面 的 表面 能 和 功 函数 计算 结果 , 以 及 与 文献 报道 的 计算 结果 和 实验 数据 "的 比较 


Table 2 Surface energy and work function of (0001) surface, and comparisons of calculated and experimental results" 


Metal Surface energy / (J-m ) Work function /eV 
This work Calculation Experiment This work Calculation Experiment 
Be 1.76 1 83 1:020Y “1.63 2700 5.23 5.62™ 4.980" 
Mg 0.57 0.64®, 0.79"" 0.78"™, 0.76®" 3.71 3.88™ 3.66®! 
Ti 1.96 1.46®, 1.82°Y 1.99"”, 2.10™ 4.25 4.59™ 4.332" 
Zr 1:59 1.73P 1.53® 1.91", 2.00°" 4.16 4.15™ 4.05P 
3.0 0.9 
(a) (b) 
2.8| es (0001) “和 
a 027 Ey 
E 2.6 on E oe 
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Fig.3 Calculated surface energies of Be (a) and Mg (b) with different surface indices and slab thicknesses 
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与 功 函 数 的 高 低 顺序 相反 , 由 低 到 高 为 : X0001)< 
M0221)<X1010)<X1231)<X1230)<X(O111). 
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同一 主 族 , 可 以 认为 原子 序数 小 的 一 方 其 功 函 数 和 
表面 能 会 更 高 一 些 , 这 与 Skriver 等 中 和 Wang 等 中 的 


功 函数 反映 了 表面 得 失 电 子 的 难 易 程 度 ， 
(0001) 表 面 功 函数 最 高 , 且 其 表面 能 最 低 , 说 明 其 最 
难 失 去 电子 , 表面 稳定 性 最 高 ; (0111) 表 面 功 函数 


研究 成 果 一 致 . 对 于 同族 元 素 , 原子 序数 越 大 , 原子 
活泼 性 越 高 , 同一 表面 功 函 数 就 会 低 一 些 , 但 其 表 
面 能 会 相对 原子 序数 低 的 元 素 小 一 些 . 这 是 因为 功 


在 几 个 表面 取向 中 最 低 , 而 表面 能 又 最 高 , 反映 了 
这 一 取向 的 表面 稳定 性 很 差 , 在 实际 情况 中 很 难 稳 
定 存在 . 对 于 金属 晶体 , 在 承受 切割 、 形 变 等 操作 
时 , 其 暴露 在 金属 最 外 层 的 唱 面 是 很 难 控制 的 , 而 
不 同 表面 的 功 函数 表面 能 差异 很 大 , 稳定 性 差异 很 
大 , 会 影响 材料 的 使 用 性 能 . 

表面 能 反映 了 表面 的 能 量 高 低 , 直接 关系 到 材 
料 表面 的 稳定 性 和 发 生 重 构 的 可 能 性 等. 计算 得 到 
的 几 种 金属 不 同 表面 的 表面 能 变化 较 大 , 且 表 面 能 
与 其 功 函数 之 间 还 存在 着 逆向 关系 . 对 于 材料 的 府 
蚀 而 言 , 应 当 尽量 避免 表 面 能 高 而 功 函 数 低 的 表面 
暴露 在 金属 外 部 , 对 于 存在 某 一 指数 表面 的 功 函数 
很 低 的 金属 也 要 做 相应 的 防腐 措施 

从 图 4 还 可 以 看 出 , Be 的 功 函 数 要 比 同一 表 二 
取向 Mg 的 功 函数 高 一 些 ; Be 的 各 表面 的 表面 能 也 
要 比 相 应 表面 Mg 的 表面 能 高 一 些 . 鉴于 两 者 处 于 


函数 反映 了 原子 得 失 电子 的 能 力 , 同族 元 素 中 原子 
序数 越 小 , 原子 半径 越 小 , 则 原子 核对 核 外 电子 的 
控制 越 强 , 电子 越 难 从 金属 表面 逃逸 , 功 函 数 也 就 
越 大 . 而 表面 能 反映 了 原子 表面 能 量 的 高 低 , 原子 
序数 越 小 , 晶 胞 的 参数 也 就 越 小 , 表面 原子 受到 周 
围 其 他 原子 的 作用 越 强烈 , 则 表面 能 量 容易 变 高 . 
3.2 Ti 和 Zr 的 计算 结果 与 讨论 

对 于 Ti 和 Zr2 种 金属 进行 了 类 似 于 Be 和 Mg 
的 表面 相关 计算 . 结果 如 图 5 所 示 . 可 以 看 出 , Ti 和 
Zr 的 功 函 数 和 表面 能 的 计算 结果 显示 的 规律 与 Be 
和 Mg 很 相似 . 对 于 Ti 而 言 , (0001) 表 面 的 功 函 数 最 
高 , 为 4.25 eV, (0111) 表 面 的 功 函 数 最 低 , 为 3.22 eV, 
不 同 表面 功 函数 差异 达 1.03 eV; (0111) 表 面 的 表面 
能 最 高 , 为 2.51 Jjm2, (0001) 表 面 的 表面 能 最 低 , 为 
1.96 J/m, 表面 能 差异 达 0.55 J/mr. 对 于 Zr 而 言 , (0001) 
表面 的 功 函 数 最 高 , 为 4.17 eV, (0111) 表 面 的 功 函 数 最 
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低 , 为 3.20 eV, 不 同 表 面 功 函数 差异 可 达 0.97 eV; 
(0111) 表 面 的 表面 能 最 高 , 为 2.09 Jnm2, (0001) 表 面 的 
表面 能 最 低 , 为 1.60 J/m’, 差异 达 0.49 J/m. 

同样 , 两 者 各 表面 功 函数 也 有 着 规律 的 排 
列 , 由 高 到 低 近 似 遵 从 : 2(0001)> 23(0221)>@(1010)> 
G1231)>G(1230)>G(011D). 表面 能 的 排列 规律 为 : 
M0001)<X0221)<X1010)<X1230)<X1231)<X0111). 
除了 X1230) 和 X1231) 表 面 顺序 变化 外 , 其 它 表面 的 表 
面 能 大 小 顺序 与 Be 和 Mg 的 表面 能 顺序 完全 相同 . 

计算 结果 还 表明 , Ti 的 功 函 数 要 比 Zr 的 大 一 
些 , 这 与 Be 和 Meg 的 功 函 数 关 系 一 样 . 表面 能 
存在 与 Be 和 Mg 相似 的 关系 . 
3.3 讨论 

利用 第 一 原理 方法 计算 研究 金属 材料 表面 的 功 
函数 和 表面 能 , 能 够 比较 容易 地 对 一 系列 复杂 表面 
进行 系统 研究 , 其 计算 结果 与 实验 值 差 别 不 大 , 完全 
在 可 以 接受 的 误差 范围 内 , 甚至 对 于 一 些 简 单 的 唱 
面 取向 , 其 计算 结果 与 实验 值 很 接近 , 可 以 很 好 地 用 
来 探究 材料 表面 的 相关 理化 性 质 , 并 用 于 指导 实践 . 

从 4 种 金属 的 计算 结果 可 知 , 对 于 六 方 唱 系 
的 Be Mg, Ti 和 Zr 而 言 , 其 功 函 数 遵从 着 0001)> 
@(0221)>@(1010)>@(1231)>0B(1230)>9X0111) 的 排 
列 规律 , 而 表面 能 也 基本 遵从 着 X 0001 )<X0221)< 
X1010)<X1231)<X1230)<X0111) 的 大 小 排列 规律 . 
这 种 排列 规律 对 于 4 种 金属 基本 都 符合 . (0001) 和 
(0221) 表 面 的 功 函数 较 高 但 表面 能 较 低 , 而 (1230) 
和 (0111) 表 面 的 功 函 数 较 低 但 表面 能 较 高 . 他 们 之 
间 的 差 值 分 别 可 以 达到 1.0 eV 和 1.0 J/rY 左 右 , 这 
样 大 的 差异 对 于 材料 的 稳定 性 是 很 不 利 的 . 依照 
计算 得 到 的 规律 , 可 以 很 容易 推测 哪些 表面 取向 
更 不 利 , 面 对 复杂 环境 时 更 容易 重 构 、 分 解 等 , 这 
可 以 在 一 定 程度 上 改善 人 们 对 于 材料 表面 的 认识 . 

计算 显示 , Be 的 功 函 数 和 表面 能 要 比 Mg 高 一 
些 ; 下 的 功 函 数 和 表面 能 要 比 Zr 高 一 些 . 这 可 以 深 
化 人 们 对 于 元 素 周期 表 的 认识 . 不 单 是 一 些 简 单 的 
理化 性 质 会 遵从 元 素 周 期 表 的 排列 , 甚至 对 于 一 些 
不 常见 的 唱 面 取向 的 功 函数 和 表面 能 , 其 也 可 能 会 
遵从 元 素 周 期 表 的 规律 . 
4 结论 

(1) 对 于 4 种 六 方 晶 系 的 纯 金属 Be, Mg, Ti 和 Zr 的 
(0001), (0221), (1010), (1231), (1230) 和 (0111) 表 面 的 
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功 函 数 和 表面 能 的 计算 研究 表明 , 不 同 表面 的 功 函 
数 基 本 遵从 大 小 顺序 : 0001)>G(0221)>GC1010)> 
GO(1231)>GO(1230)>GCCO111D); 其 表面 能 遵从 规 得 
XM0001)<x(0221)<X1010)<X(1231D)<X1230)<xXO11 

(2) 同族 金属 中 , Be 各 表面 的 功 函 数 和 表 卫 
要 比 相 应 表面 Mg 的 大 一 些 ; 五 相应 表面 的 功 函 
和 表面 能 的 也 要 比 立 的 大 一 些 . 

(3) 所 研究 的 6 种 金属 表面 中 , 不 同 表面 取向 2 
间 的 功 函 数 和 表面 能 差异 明显 , 有 些 表面 之 间 功 函 
数 差异 达到 1.0 eV, 有 些 表 面 之 间 的 表面 能 差异 超 
过 了 0.5 Jim2. 其 中 , (0001) 表 面 最 稳定 , 而 ( 0111 ) 最 
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